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Abstract 

We clarify the path integral calculation, recently suggested by Golubev and 
Zaikin [|^J2|, and show, contrary to their claim, that quasiparticles become 
fully coherent quantum particles in the T ^ limit. The important physical 
point is the inclusion of the recoil of the quasiparticle when interacting with 
fluctuations in the rest of the Fermi gas. 
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I. INTRODUCTION 



Golubev and Zaikin (GZ) [^,0 has recently made a rather dramatic proposal, that the 
quasiparticles of an electron gas, will loose their quantum coherence in a finite — not terribly 
long — time, even at temperatures going to the absolute zero, due to the interaction with the 
rest of the gas. Obviously such a fact would have far reaching consequences. One would have 
to rethink a number of fundamental properties of electron physics, in particular the theory 
of localization. And it is indeed to good to be true. In a very careful and pedagogical work 
Aleiner, Altshuler and Gershenson (AAS) have redone the theory of weak localization, 
this time taking electron-electron interactions into account. And they show convincingly, 
that electron-electron interaction does have an effect at T = 0, but it is not to decohere the 
quasiparticles, but rather to give extra scattering off the static Friedel oscillation pattern of 
the screening cloud associated with the impurities. 

Still the GZ work is not without merits. It wants to formulate the many body transport 
theory problem using Feynman path integrals. This has several advantages. The physical 
picture emerging is much more readily understandable. The physics of weak localization 
itself was not widely known and accepted until the works of Altshuler et al. 0, Bergmann 
0] and Chakravarty and Schmid |^ , where a real space description in terms of electron paths 
was used. The generalization of such an approach to interacting many body systems is still 
not fully understood — the GZ work is certainly an example of this. It therefore is of interest 
to find out, where exactly the GZ calculation goes wrong. The AAS calculation — using 
standard methods of many body theory — does show that it is wrong, but it does not really 
show how the procedure can be modified. The only hints that AAS offers is some vague 
statements about problems in semiclassical calculations when classical orbits intersect. 

In this paper we want to show that the error in the GZ calculation is not so much related 
to the selfcrossing orbits relevant for weak localization, but rather stems from the fact that 
the recoil of the particle, who's path is being followed, in the real and virtual scattering 
processes with the rest of the gas, is not properly taken into account p. In section 2 we 
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show how to set up a more reahstic path integral calculation in the Caldeira-Leggett spirit 
also used by Golubev and Zaikin, and in the final section we show how scattering processes 
are properly dealt with, giving the result, that in the T = limit, the dephasing time for 
quasiparticles at the Fermi level will tend to oo. 



II. PATH INTEGRALS OF AN INTERACTING FERMI GAS 

Our starting point is the same as that of Golubev and Zaikin. We want to consider an 
interacting electron gas, i.e. to consider a many body problem with the following Hamilto- 
nian: 

H — Ho + Hint, (1) 



where 



Ho^ jdr ^lir) 



Mr): (2) 



and 



Hint ^Ijdrjdr' iPl{r)i^Ur')v{r - r')^.'(^')^a(r) (3) 

where U{r) is an external potential which may or may not be random, n is the chemical 
potential, and v{r — r') is the interaction potential which we shall take to be the Coulomb 
potential — others may be used, that will not affect the points we want to make in this 
paper. 

The goal is to calculate the many body density matrix, and from that get physical 
quantities like currents and densities. We shall follow GS and use a functional integral 
formulation, introduce a Hubbard-Stratonovic field Vi{r,t) {i being a Keldysh index), to 
decouple the electrons. The effective action for the fluctuating V-flelds is calculated by GS 
to second order. Since the approximation of leaving out higher order terms is not being 
questioned, we shall do the same, and hence work with 
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tS[Vi,V2]=tJ -j—^ i^V {-uj,-k) V+{uj,k) 

+^ V-{-uJ, -k) i^r^coth i^—j V-{u, k)j . (4) 

Here V'^ — {Vi + V2)/2 and V~ — Vi — V2. It is important to notice the inverse of the 
permitivity is a retarded function so for instance is the average {V~ {—k,t')V^ {k,t")) zero 
for t' > t" . The physical density matrix p{t) is obtained by averaging the density matrix in 
the presence of specific field history pv{t) over V 



Pv{t) is given by 



with the effective electron action 

S^, tP, V;t]^ 1^ dt' dr r)^t>^P{t\ vr) - e^{t', r)dt>tP{t\ vr)] - //q) (7) 

The Keldysh contour Ct runs from —00 up to time t and then back to —00. By direct 
differentiation a differential equation for pv(t) can be derived and the linear response to an 
applied electrical potential eVx is then given by 

.95^ = H,5pv{t) - dpv{t)H2 - [eV,, pv] (8) 

The solution to this equation is (GZ Eq. (41)) 

6py{t) = I f'dt'U,{t,t')[eV,,pv{t')]U2{t',t), (9) 
Jo 

where the time evolution operators are given by 

f*2 



Ui,2{ti,t2) = Texp 



t dt'Hi,2{t') . (10) 
Jti J 



and the effective Hamiltonians Hi 2 are functionals of the fluctuating Hubbard-Stratonovic 
fields and are given by 



H^ = Ho-eV+--il-2po)eV- 

H2 = Ho- eV+ + leV~{l - 2po), (11) 
where po is the equihbrium one particle density matrix 

This whole theoretical setup is very close to the standard Caldeira-Leggett description of 
a single quantum particle in a dynamic environment, as also emphasized by Golubev and 
Zaikin in a more recent paper where they reply to some of the critizism raised. The 
environment is responsible for both a kind of "dynamic Hartree" contribution —eV~^ to 
the effective Hamiltion and a non-Hermitian contribution — 1(1 — 2pQ)eV~ (to Hi, and a 
similar to H2). In contrast to the standard one-particle Caldeira-Leggett, the non-Hermitian 
contribution is very sensitive to the energy of the quasiparticle, through the factor 1 — 2po, 
which essentially is a smoothed out sign function. This of course has to do with the Pauli 
principle playing an important role in the many body problem. 

At this point GS represent the one-particle time-evolution operators Ui and U2 as stan- 
dard real space path integrals, i.e. the density matrix (^ is an integral over pairs of paths. 
Quantum interferense effects, like the ones responsible for weak localization, only occur if 
the two paths in the pairs are allowed to differ substantially. Such pairs of different paths 
are suppressed due to the coupling to the V^s in the effective Hamiltonians. If the pairs 
differ, they will experience uncorrelated V"^ so, in order to study the strength of dephasing 
one only need to study one of the time-evolution operators Ui^2 and see how fast it decays. 

So we want to make a path integral representation of Ui{t). As mentioned, the non- 
Hermitian part of the Hamiltonian is strongly energy dependent, so a path integral, where 
one has a lax attitude to the exact energy of the paths, as do GS when they in their 
paper make the approximation, that the relevant paths are semiclassical orbits (straight 
lines between impurity collisions) with constant energies. On the scale of the Fermi energy, 
ep this may be a good approximation, but a change of energy of order fc^T in a electro n- 
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electron collision event may have a dramatic effect in the contribution from non-Hermitian 
term, which exactly varies on this much smaller energy scale. Mathematically the ordering 



of po and V~ in (|TT]) is very important, but this aspect is missed in a conventional path 

integral representation with classical orbits. 

To repair the calculation we on the one hand want to follow the eletrons in real space, 

since this is at the heart of e.g. weak localization, and on the other hand keep track of the 

energy and momentum to a precision much smaller than the Fermi energy. We therefore 

represent the one particle time evolution operator Ui{t) as a path integral where we instead 

of using position eigenstates use wavepackets \x,p): 

/ , V f {r-xY . \ 
{r\x, p) oc exp I — — + ip ■ r \ , 

real space wavefunction 

{q\x,p) oc exp I ^ + i[p - q) ■ x \ , 

momentum space wavefunction, (13) 

where the width a will be chosen appropriately, e.g. 



Af/, (14) 

where Ai? is the Fermi wavelength and / is the eleastic mean free path. In this way we can 
locate the electron so well in space that it makes sense to talk about orbits that travel from 
one impurity to another. On the other hand the momentum (and hence the energy) is very 
well determined on the scale of the Fermi momentum and we can keep track of small changes 
in energy. 

In this basis of wavepackets free propagation of an electron is given by 

(p',x'|e-*^°*|a;>P) 

{x' — X — vty {jp — p'Ya^ 



exp 



exp 



2ma 



4(a2 

/ IrnvH -k-{x + x')/2 + (^)' ^ + p ■ x - p' ■ x')\ 



—I- 



\ 



, 2 

1 ■ ' 



t 



2mo? 



(15) 
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where 

= k^p-p. (16) 

In the following we shall use the simple approximation 

{p', x'\e-'"°^\x,p) ^ d{x' -x-vt) S{p - p') exp {-iept) , (17) 

which of course is nothing but semiclassical propagation. 

The matrixelements for the coupling to the rest of the electron gas is given by 

{p',x'\eV+{r,t)\x,p) 

= J dkeV+ik, t) exp (-^^^^ - + - z(p + fc - p') . {x + x')/2^ 

fti 6{x - x') J dk 6{p + k-p) V+{k, t) (18) 

and 

{p',x'\^{l-2po)eV-ir,t)\x,p) 

p' — p — kYa^ {x — x'Y , . x-\-x' 2/ .{x' — x)' 



J f J, / [p -p - K) [x-x'y , . x^x' 2/ 
dq J dk exp ^ ^^2 - ^(P - P-k)- a'U-i 



2a? 



^(1 -2po(e p/+p+fc ^P)ey {k,t) 



5{x - x') I dk 6{p' -p-k)^{l- 2po{ep+k))eV-{k, t). 



We now carry out the standard path integral construction of slicing the time into N 
infinitesimal pieces 5t, and insert an (over-)complete set of wavepackets between the Ui{5t) 
factors: 

7 = 1 W(p/,a;;|f/i(t)|x„p,)e^^[^i'^^l 

= /^^n / dXnj dpn{Pn+l,Xn+l\e-'''^'-')''\Xn,Pn)e''^^"'''^. (20) 

In this expression we want to integrate out both momenta and the fluctuating V fields, 
leaving us with a real space path integral. To obtain the effective action for that integral, 



we use a variation of Pauli's standard A trick. First the interaction {H' = H — Hq) is 
multiplied by a factor A. Then we consider a functional g{X, x) of the real space orbits: 

g{\x) = I VV\{ I f/p.(p„+i,a^„+i|e-^^'(*")^*|a;„,p„)e^^[^-^^l (21) 

This functional satisfy a simple differential equation in A: 



d\ Jo {{U\t,0,x))) 

where we have taken the N ^ oo limit, introduced the notation 

{{■)) = ! VV J dp .e'^^^''^'\ (23) 

and the time evolution function 

U^{t2,h,x) = n {Pn+l,Xn+l\e-''' \Xn, Pn) , Ui = U/tN. (24) 



n=ni 



The equation (0) is easily solved, and we end up with the following expression for the path 



integral (pOD 



Px^(0,x)e^^-//(^\ (25) 



where the effective action is 



9 (x)- f'rIX f\. A{U\t,t',x)H'{t',x)U\t'Ax))) 

Like GZ we are going to evaluate this action to second order in the interaction, i.e. expanding 
the time evolution operators to only first order in A, since Seff contains an explicit interaction 
term. The final result becomes 

Q (^^ ■ f f ,. A{U%t,t\x)H'it')U%t',t'',x)H'it'')U%t'',0,x))) 

S^,j{x) =rj^dt dt WWM)) • ^ ^ 

This action is a functional of the real space path being considered. For simplicity we shall 
evaluate it for straight line path with a constant momentum p very close to the Fermi 
momentum. Using the semiclassical propagator ( P^ ) and the matrix elements (|TB|) and (|TUp 
this action is readily evaluated to 
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Seffix) = f^dt' dt" I dk (^{V+{k,t')V+{-k,t")) + ^(1 - 2po{ep^k)){V^{k,t')V'{-k,t")) 



g-i{ep_fe-ep)(t'-t") 

duj 
2^ 



dk 



(^{V+{k,u;)V+{-k,-u;)) + ^{l-2po{ep_k)){V+{k,u;)V-{-k,-u;))'^ 



f{-uj + v-k,t). (28) 

Here we have used that {V^{—k, t')V^{k, t")) is an advanced function and hence zero in our 
context. The function f{Q,t) carries the time dependence and is given by 
N /■* , / /■*' , // io(f' i"] . /I sin fit \ 1 — cos fit 

At this point we pause and make a detailed comparison to the GZ calculation. The last 
term of (pSD is equivalent to the first term in GZ eq. (54), taking into account, that we only 
calculate the contribution to the action from one path. If the GZ function R is written out, 
one obtains in the GZ version (using r{ti) — r(t2) = v(ti — t2)) 

and in our version, with (y^{—k, —uj)V^{k,uj)) written out 

where we now use the GZ notation n{p) = po{tp). The difference is, that in our formula the 
Pauli factor depends on the recoiled momentum p+k whereas in GZ this factor is depending 
on p\ This we believe is precisely where the GZ calculation goes wrong. 

We shall now continue with our calculation, and show that it gives the expected result, 
i.e. no dephasing for electrons at the Fermi level at T = 0. 

First, for /(fi,t) in the long time limit we have 

/(fi,t)= fdt' r'citv^(*'-*") 

Jo Jo 

/ 1 sin fit \ 1 — cos fit 

U""fi^j +^ fi^t 

~ it ^— — ?7r(5(fi)^ large t 
it 



fi + i?7 



rj a positive infinitesimal. (32) 



Since {V~^{—k, —uj)V~{k, u)) is the Fourier transform of a retarded function, it is analytical 
in the lower complex uj plane, and the u integral in (pSf ) is performed by closing the contour 
in the lower half plane and picking up a contribution from the pole at u; = —k ■ v — irj. 

Since our main concern is with the loss of coherence, we shall only evaluate the imaginary 
part of Seff- After some simple manipulations we get our final and main result in the long 
time limit 

-'-K^)) '''' 

In the zero temperature limit, which is our main concern this reduces to 

ImS.fjix) = 1 1 dk^-^^^^^^^i4 - ■ ^)^)sign(^ ■ k). (34) 

This is by no means surprising. We actually have recovered the standard result from standard 
many body theory. The coefficient to t is nothing but the imaginary part of the self energy 
of a particle with momentum p. This does go to zero when the momentum approaches 
the Fermi momentum, resulting in infintely long lifetimes — and coherence times — for a 
quasiparticle at the Fermi surface! 

III. CONCLUSION 

The claim by Golubev and Zaikin, that quasiparticles at the Fermi level will decohere 
due to the interaction with the zero-point fluctuations of the modes of the rest of the Fermi 
gas has been shown to be due to an error when dealing with the real space path integral 
representation of quasiparticle propagation. If one does not allow for tiny changes in the 
momentum and energy — not significant for the overall path, who's shape is dominated by 
the fast velocity, vp — then the erroneous result occur. The cure is to use wavepackets, with 
reasonably well defined space and momentum coordinates. The space part of the resulting 
path integral takes care of the overall propagation of the particle, whereas the momentum 
coordinate is important when calculating the interaction with low energy modes of the rest 
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of the system. This interaction is very much dominated by the Pauh exclusion principle, 
which is extremely sensitive to minute changes (on the scale of ep) in the energy. The latter 
effects are of course much more precisely described in standard many body theory, which 
usually uses a basis of momentum eigenstates. 

From our point of view the story is now over. We have shown that quasiparticles can 
indeed stay coherent and e.g. give an interference pattern in a double slit experiment. In 
such an experiment the paths going through one or the other slit are very different, making 
our considerations with only one path relevant. GZ are studying weak localization, where the 
interfering paths are time reversed paths, which actually visit the same spatial points (but 
at different times), giving an unlikely loophole for such paths to decohere. In an appendix 
we show that also these time reversed paths interfere! 
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APPENDIX: 

For the forward in time propagator {/i we found the effective action due to the interaction 

Seffix) = l^dt'J^ dt" J dk (^{V+{k,t')V^i~k,t")) + ^(1 - 2po{ep.k)){V+{k,t')V-i-k,t"))^ 

= jdk J^i^{V+ik,uj)V+{-k,-uj)) + ^il-2poie,^k)){V+ik,uj)V-i-k,-Lu))^ 
f{-uj + v-k,t). (Al) 

Likewise we find for U2 the effective action 
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suffix) = J^dt'J^ dt"J dk (^{V+{k,t')V^{-k,t")) - ^{l-2po{e,+k)){V+{k,t')V-{-k,t")) 

-Jdkj^ (^{V^{k,u)V^-k, -co)) - 1(1 - 2poie,+k)){V^k,co)V'i-k, -a;))) 
f{-u + v-k,t). (A2) 

The pairs of paths that contribute to the weak locahzation effect are the ones where forward 
1 and backward in time 2 paths are each others time reversed partners. In this case it is 
also possible to get cross contributions where one V is on the forward path while the other 
is on the backward in time path. These contributions correspond in GZ to the two last 
terms in (54) and (55) and after focusing on the weak localization (68) and (69). Like GZ 
we find that the cross contributions disappear after averaging over the diffusive paths. The 
remaining term in GZ (69) originates from the average in the above actions and is 

identical to what GZ find in (71). In GZ the first term in (68) vanishes. We find in analogy 
with the above single line calculation an extra contribution proportional to tanh (^|^)- 
To see this explicitly let us consider a particular straight line segment of the path and let 
us consider the case where the interaction also takes place within this straight line segment. 
If the velocity along the forward path 1 is v is the velocity along the time reversed path 2 
—V. The combined contribution from the second term in the effective actions is 

jdkj ^{po{ep+k)-Po{ep-k)){V+{k,u;)V-{-k,-u;))f{-u; + vk,t). (A3) 

We have used {V'^{k,u;)V''{—k, —u;)) is even in k and e_p_fc = e^+fc 
Compared to GZ (68) we see that the only difference is that instead of a factor (poi^p) — 
Po(ep)) = we get 

tanh - tanh (^^) 

(po(6p+.) - po(6p-.)) = y^h^L±_ — yjtB7i2^ (A4) 

so in GZ the coth (2^^) should always be replaced with 



coth 



UJ 



tanh C^^-^) , (A5) 



2kBTj V 2kBT 



which vanshes in the T — > and — > limit. 
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